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Internal colour gradients for E/S0 galaxies in Abell 2218 
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We determine colour gradients of —0.15 ±0.08 magnitudes per decade in radius in 
F450W-F606W and -0.07±0.06 magnitudes per decade in radius in F606W-F814W 
for a sample of 22 E/S0 galaxies in Abell 2218. These gradients are consistent with the 
existence of a mild (~ —0.3 dex per decade in radius) gradient in metal abundance, 
(cf. previous work at lower and higher redshift for field and cluster galaxies). The 
size of the observed gradients is found to be independent of luminosity over a range 
spanning M* — 1 to M* + 1.5 and also to be independent of morphological type. These 
results suggest a fundamental similarity in the distributions of stellar populations in 
ellipticals and the bulges of lenticular galaxies. These results are not consistent with 
simple models of either monolithic collapse or hierarchical mergers. 
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1 INTRODUCTION 

In Hubble's (1936) classification scheme for galaxies, lentic- 
ular (SO) galaxies first appear as a transitional class between 
the ellipticals and the two branches (barred and unbarred) 
of the spirals, and, despite much evidence to the contrary, 
it has often been assumed that the Hubble 'fork' describes 
an evolutionary scheme. For instance, this assumption in- 
forms the popular idea that collisions of proto-disks in the 
early universe lead to the formation of elliptical galaxies 
or that quenching of star formation in field spirals as they 
fall into dense regions leads, by a variety of mechanisms, 
to their transformation into SO galaxies. This process, first 
proposed by Larson, Tinsley & Caldwell (1980), has been 
invoked to account for the morphology- density relation of 
Dressier (1980) and the apparent blueing of the galaxy pop- 
ulation in distant clusters (Butcher & Oemler 1984), and 
has stimulated considerable theoretical activity (see Pimb- 
blet 2003 for a review). 

Such processes modify the star formation history of the 
disk by removing the supplies of gas necessary to fuel further 
star formation. A sensitive test of this mechanism would 
therefore involve a comparison of the stellar populations of 
bulges and disks for elliptical, SO and spiral galaxies. While 
the stellar content of galaxies cannot be resolved except for 
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some of the nearest objects, the radial distributions of their 
stellar populations may be examined by means of colour 
gradients. 

It has long been known that early-type galaxies exhibit 
negative colour gradients, in the sense of being bluer at 
larger radii (Franx, Illingworth & Heckman 1989, Peletier 
et al. 1990, Goudfrooij et al. 1994, Tamura et al. 1999, Saglia 
et al. 2000, Tamura & Ohta 2000, Hinkley & Im 2001, Idiart, 
Michard & de Freitas Pacheco 2002, Tamura & Ohta 2003). 
These gradients have been conventionally interpreted as rep- 
resenting gradients in mean metal abundance, as most local 
ellipticals also exhibit such gradients in the strength of metal 
absorption lines (Baum, Thomsen & Morgan 1986, Carollo, 
Danziger & Buson 1993, Davies, Sadler & Peletier 1993). 
However, this intepretation is vitiated by the existence of 
a strong degeneracy between age and metallicity in broad- 
band colours as well as in most spectrophotometric indices 
(Worthey 1996). Nevertheless, colour gradients have con- 
siderable potential for constraining the history of star for- 
mation in distant galaxies, and for providing clues to their 
formation and evolution history. 

One approach to breaking the age-metallicity degener- 
acy is to consider the evolution of colour gradients. This ex- 
periment has been attempted in a few distant clusters with 
rich populations of early- type galaxies. Saglia et al. (2000), 
Tamura & Ohta (2000) and La Barbera et al. (2003) have 
observed galaxies in clusters at z ~ 0.4 and studied the evo- 
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lution of their colour gradients, concluding that the most 
likely explanation consists of a mild metal abundance gra- 
dient within a passively evolving stellar population formed 
at very high redshift. 

However, because of the relatively limited number of 
galaxies and relatively poor resolution at these redshifts, 
even with Hubble Space Telescope (HST) imaging, Saglia 
et al. (2000), Tamura & Ohta (2000) and La Barbera et al. 
(2003) have not examined the dependence of colour gradi- 
ents on galaxy luminosities and morphology and whether 
the disk and bulge components of galaxies have different 
gradients, and therefore different stellar populations. It is 
therefore useful to consider a somewhat lower-redshift tar- 
get, where the galaxy population can be probed over a wider 
range of luminosities and where the higher spatial resolu- 
tion (because of the lower distance) and smaller cosmologi- 
cal dimming allow more detailed studies of colour gradient 
evolution than has been possible (cf. Tamura & Ohta 2003 
for a similar approach to the z = 0.03 cluster Abell 2199). 

In this paper we discuss a study of colour gradients in 
the z = 0.18 cluster Abell 2218 from deep archival HST 
multicolour data. We describe the data and the analysis in 
the following section and present the main results. We dis- 
cuss these results in the light of models of galaxy forma- 
tion and, particularly, disk evolution. We adopt the concor- 
dance cosmological model with Q,m = 0.3, Qa = 0.7 and 
H = TOkms^Mpc- 1 . 



2 DATA ANALYSIS 

The data consist of three archival images of the core of Abell 
2218 taken by the HST with the Wide Field and Planetary 
Camera 2 (WFPC2) through filters F450W, F606W and 
F814W (PID: 8500; PI: Fruchter). For simplicity, we will 
refer to these bands as B' , V' and I' in the remainder of 
this paper (we use prime symbols to distinguish the bands 
from the actual Johnson bands and to avoid confusion in 
the subsequent discussion). The total exposure times were 
12000s, 10000s and 12000s respectively, with individual ex- 
posures being 1000s each. These images were retrieved as 
fully processed associations from the HST archive ( Micol, 
Bristow & Pirenne 1997). A more detailed description of this 
dataset is given by Smail et al. (2001). 

We have selected a sample of 24 galaxies morphologi- 
cally classified as E or SO and spectroscopically confirmed 
to be cluster members (this information is taken from Smail 
ct al. 2001). First, the mean sky was removed from the im- 
ages using the IRAF task sky. For each galaxy we then 
derived surface brightness profiles using the IRAF task 
ellipse (Jedrzejewski 1987, Busko 1996), using 0.1" (1- 
pixel) sampling and 3<r clipping to remove overlapping stars 
and galaxies. As our main interest in these data is to de- 
rive colour gradients, we forced the isophotes in V' and 
I' to be evaluated at the same radius (measured along the 
semi-major axis) as in the B' image (which is the one with 
the lowest signal to noise) by using the 'inellip' option in 
ellipse. The isophote centres were left free to be recentered 
for each image, to avoid the effect of minor misalignments 
between the HST images in the different bands. 

In some cases, galaxies lie in the envelope of other bright 



galaxies or overlap with nearby galaxies (especially the two 
giant ellipticals but also in a few other cases). In these cases, 
we have modelled the 'contaminating' object using the IRAF 
task bmodel and removed it from the image. We then com- 
puted new isophotes for the galaxies. In most cases, this 
appears to make no difference, at least to the level of pro- 
ducing colour gradients different by more than the statistical 
and systematic errors. However, galaxies # 292 and #298 
appear to be significantly affected by the envelopes (and, in 
the case of # 292 proximity to the edge of the detector) of 
their neighbouring galaxies. Because of this we decided to 
exclude these two galaxies from our analysis. In other cases, 
such as # 205 and # 634, modelling shows that the profiles 
are not being affected by proximity to a bright neighbour 
and that ellipse has appropriately removed the contami- 
nating object from the profile.. 

For each data point, two sources of error contribute; us- 
ing the terminology of Saglia et al. (2000) these are: a sta- 
tistical error, which is the uncertainty, returned by ellipse 
in measuring the mean flux along each isophote and a sys- 
tematic error, reflecting the uncertainty in measuring and 
subtracting the mean sky level. We estimated this latter 
component in the following way: for each of the individual 
1000s exposures in each band, we calculate the mean sky us- 
ing sky. This yields 12, 10 and 12 sky estimates in B' ,V' and 
/'. The error in the mean sky is then computed by jackknifc 
resampling of these sky estimates. 

We calibrate the data on to the Vega system, using 
the zeropoints provided by Holtzman et al. (1995). At each 
isophote semi-major radius we then compute the resulting 
colour as a function of radius. 

We measured the point spread function using non- 
saturated stars in our images. There are, unfortunately, only 
a few such stars but the full width at half maximum appears 
to be a consistent 0.15" in all bands, suggesting that differ- 
ences in the point spread function do not affect the deriva- 
tion of colour gradients. Inspection of the surface brightness 
profiles in Figure 1, shows some deviation from the profile 
expected for a de Vaucouleurs profile in the inner two or 
three points. For this reason, we decide to fit to the points 
with r > 0.3" in order to avoid any residual effects due to 
small differences in the point spread function. 

The colour gradients are then derived by fitting the data 
(colour vs. logarithm of the radius measured along the semi- 
major axis) with a weighted linear least squares regression, 
where the weights are the statistical errors of the colours 
(added in quadrature). The fits are carried out to either 
the radius at which the quadratic sum of statistical and 
systematic errors exceeds 10% or to the radius at which 
ellipse stopped integrating the surface brightness profile 
(even if the errors are smaller than 10%). We also estimate 
the colour gradient only over the range defined by our worst 
colour (£>' — V), for consistency. The statistical error on 
the fit slope and intercept is derived from the error in the 
fit to the data, while the systematic error is determined by 
carrying out a new fit to the data points, after adding and 
subtracting the ±la error in the sky level. 

Figure 1 shows the surface brightness profiles in all three 
filters and the colour gradients in B' — V' and V' — I', plot- 
ted against logr, for all 22 galaxies. The surface brightness 
profiles and colour profiles are plotted over the range r = 
to the largest radius used in fitting the colour gradients. 
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Figure 1. Surface brightness profiles and colour gradients for galaxies in A2218 (as identified in the figure). In each figure panel (one per 
galaxy), the top panel shows the surface brightness profiles in B' (circles), V (squares) and /' (triangles) plotted against logr (measured 
in arcscconds along the semi- major axis); the middle panel shows the colour distribution as a function of logr in B' — V' (data points 
are the filled circles, error bars are the statistical error in the flux, the dashed lines show the effect of a ±lcr error in subtracting the sky 
level, and the thick dashed line shows the best fitting straight line); the bottom panel is the same as the middle panel, but for V — I' 
and using squares. The open symbols in all panels show the surface brightness values and colours for points interior to r = 0.3" which 
were not used in the fit. Note that we refer to the B' , V and /' colours by their proper HST names, to avoid confusion for those who 
may wish to make use of these data (the actual profiles are available on request in machine-readable format). 
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Figure 1. (continued) 



As for points interior to 0.3", which are not used in fitting, 
some galaxies show colours consistent with the extrapolated 
fits, while in some other cases it is possible that the galaxies 
contain 'cores' with flatter colours. However, this is difficult 
to ascertain without data at better resolution. 

The major contribution to the error budget comes, not 
unexpectedly, from the B' band, where the WFPC2 camera 
is less efficient and the galaxies are fainter (as it corresponds 
to the rest-frame U). We tabulate the results in Table 1. This 
shows, in column order, the galaxy ID, morphology and K 
band luminosity (from Smail et al. 2001), the B' — V' colour 
gradient with its statistical and systematic error, and the 
and V'—I' colour gradient with its statistical and systematic 
error (derived as described above). The mean gradients are 
A(B' - V')/Alogr = -0.14 ±0.08 for Es and -0.16 ±0.09 
for SOs and A(V - I')/A\ogr = -0.06 ± 0.04 for Es and 
-0.08 ± 0.06 for SOs. 



3 RESULTS 

In order to model the stellar populations of the galaxies and 
derive gradients in age and/or metal abundance from the 
observed colour gradients, we need to determine the cen- 
tral colours as a boundary for the models. We plot these 
colours in a colour-colour diagram in Figure 2 (where the 
colour at r — is extrapolated from the linear fits), to- 
gether with the predictions for a single-age stellar population 
with present age of 12.5 Gyr and metallicity varying between 



[Fe/H]=±0.55 and —2.5, and for a single-metallicity stellar 
population with metallicity of [Fe/H]=±0.55 and ages from 
12.5 to 2.5 Gyr, as observed through the HST filters and 
at z = 0.18. We used the latest version of the GALAXEV 
models of Bruzual & Chariot (2003) to realize these simu- 
lations adopting the Padova 1994 isochrones and Chabrier 
initial mass function, as recommended by Bruzual & Chariot 
(2003) and in the GALAXEV documentation. 

One possible problem is that, as shown in Figure 2, 
the models are not fully capable of reproducing the central 
colours. This may be due to some small discrepancy in the 
population synthesis models, or to some inaccuracies in the 
extrapolation of central colours from the r > 0.3" data, or 
to the presence of dust in the galaxy cores (e.g. Peletier 
et al. 1999). However, the reddest (oldest and more metal 
rich) models are reasonably close to the actual values of 
the colours, especially in the B' — V' colour which is most 
sensitive to stellar populations and most useful for our study, 
and we adopt this as the starting point for our modelling in 
interpreting the derived colour gradients. 

One of the aims of this paper is to investigate not only 
the evolution of the colour gradients at intermediate redshift 
and discuss them in the light of both local and higher red- 
shift data (see below) but also to consider how the gradient 
size depends on galaxy luminosity and on morphology. We 
plot the derived colour gradients vs. K magnitude (which is 
a good measure of the underlying stellar mass) in Figure 3. 
We use different symbols (circles and squares, as indicated 
in the caption) for E and SO galaxies. We find little evidence 
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Table 1. Colour gradients 



ID 


Morphology 


K 


A(B' - V')/d\ogr 


Cstat 


^sys 


A(V" -/')/<* log 


r u 3tat 


asys 


301 


cD 


13.19 


-0.059 


0.002 


0.016 


-0.006 


0.001 


0.001 


307 


E 


13.77 


-0.124 


0.004 


0.017 


-0.042 


0.002 


0.001 


280 


SBOa 


14.37 


-0.180 


0.006 


0.024 


-0.092 


0.004 


0.002 


315 


E 


14.55 


—0.247 


0.007 


0.022 


-0.064 


0.006 


0.002 


1057 


E 


14.57 


-0.026 


0.006 


0.042 


-0.056 


0.002 


0.004 


428 


SO 


14.60 


-0.076 


0.006 


0.040 


-0.056 


0.003 


0.002 


401 


E 


14.85 


-0.188 


0.005 


0.044 


-0.036 


0.002 


0.003 


633 


SO 


14.88 


-0.026 


0.001 


0.025 


-0.045 


0.006 


0.002 


4007 


SB0 


15.17 


-0.125 


0.011 


0.030 


-0.075 


0.007 


0.002 


337 


E 


15.20 


-0.073 


0.009 


0.016 


-0.072 


0.005 


0.001 


205 


SO 


15.21 


-0.361 


0.015 


0.020 


-0.237 


0.011 


0.001 


137 


SB0 


15.52 


-0.126 


0.011 


0.062 


-0.022 


0.005 


0.006 


376 


SO 


15.54 


-0.123 


0.012 


0.032 


-0.102 


0.006 


0.002 


1047 


E 


15.77 


-0.061 


0.017 


0.044 


-0.064 


0.010 


0.002 


537 


E 


16.05 


-0.105 


0.023 


0.020 


-0.113 


0.016 


0.001 


4014 


SO 


16.22 


-0.124 


0.018 


0.045 


-0.092 


0.009 


0.003 


449 


E 


16.27 


-0.235 


0.010 


0.031 


-0.034 


0.005 


0.001 


154 


E 


16.54 


-0.178 


0.024 


0.075 


-0.072 


0.011 


0.006 


612 


SOa 


16.55 


-0.230 


0.018 


0.056 


-0.024 


0.007 


0.002 


131 


E 


16.93 


-0.212 


0.019 


0.075 


-0.140 


0.008 


0.003 


599 


SO 


17.32 


-0.187 


0.024 


0.065 


-0.058 


0.011 


0.004 


4003 


E 


17.46 


-0.220 


0.028 


0.071 


-0.001 


0.012 


0.005 



that gradient size, in both colours, depends on either lumi- 
nosity or morphology. The consequences of this finding are 
examined below. 

Naturally, our conclusions prevalently concern 
spheroids and the bulge components of SO galaxies as 
our data, especially in B' where we are most sensitive 
to stellar populations, do not reach much into the disk 
dominated portion of the SO galaxies. However, in some 
cases, Figure 1 shows that we are able to sample a small 
portion of the disk. This may be the case for # 280, 
633, 205 and 137. It is apparent that the colours of these 
disks are at least qualitatively similar to those of their 
parent bulges, which would suggest that they have similar 
stellar populations. Fisher, Franx & Illingworth (1996) 
and Peletier & Balcells (1996) have also reached similar 
conclusions, although Fisher et al. (1996) also point out 
that the star formation histories may be different for at 
least some galaxies, despite their having similar ages. 

Unfortunately, the small portions of the disks we survey 
are not sufficient to determine colour gradients for the disks 
and carry out a comparison with the bulge and therefore 
reconstruct the star formation and enrichment histories of 
the two components of the galaxies. This should be possible 
with more sensitive data from the ACS, some of which are 
available publicly and will be the subject of a later paper. 



4 DISCUSSION 

We now examine the implications of our findings for models 
of galaxy formation and evolution: in particular, we consider 
simplistic models of formation by dissipational collapse and 
by hierarchical mergers as our benchmarks; these compar- 



isons should probably be intended as a stimulus for more 
comprehensive modelling rather than setting true limits on 
the accuracy with which the models represent the observa- 
tions 

We first find that the size of the colour gradients is 
broadly independent of K luminosity. Models of colour gra- 
dients in galaxies, unfortunately, do not yet reach this level 
of detail, but we naively expect that, in a monolithic collapse 
scenario, the size of the gradients will be proportional to the 
depth of the potential well, as the stellar population gradi- 
ents are induced by superwinds, so that more massive galax- 
ies have steeper colour gradients (Carlberg 1984). In gen- 
eral, monolithic collapse models yield excessively large gra- 
dients when compared to the observations (Carlberg 1984) 
but more detailed treatment of gas physics may alleviate this 
problem (e.g. Kawata 2001, Pipino & Matteucci 2004). In 
hierarchical merger models (White 1980, Kauffmann 1996, 
Cole et al. 2000) we expect that mergers will erase any 
colour gradients originally set up. Some gradients may be re- 
established by star formation at later epochs, but we might 
expect that more massive galaxies, at the top of the merg- 
ing hierarchy, will experience more mergers and have flatter 
gradients. Our results appear to be inconsistent with either 
(admittedly simplistic) scenario, as we observe no evidence 
of strong trends in gradient size with luminosity. 

Tamura & Ohta (2003) suggest that the brigher (R < 
15) and more massive galaxies in Abell 2199 show steeper 
gradients, which would be consistent with the predictions of 
monolithic collapse models, while La Barbera et al. (2004) 
favour no trend in the gradient size with luminosity, or pos- 
sibly flatter gradients for brighter galaxies, in Abell 2163B, 
at z ~ 0.2 which is in somewhat better agreement with our 
observations and with hierarchical merger models. Study of 
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Figure 2. Colour-colour plot for E (circles) and SO (squares) 
galaxies and two BC2003 model lines, as identified in the legend. 
The values presented here are the central colours extrapolated 
from the fit to the data points presented in Figure 1. Solid error 
bars represent the statistical error while the dashed error bars 
show the systematic error. 



the evolution of these trends with redshift may allow us to 
resolve this issue more decisively. 

The observations suggest that there is no significant dif- 
ference in the size of colour gradients between E and SO 
galaxies bulges. This implies that the two systems have very 
similar distributions of stellar populations. We use our sin- 
gle stellar population models to calculate that the observed 
gradients are consistent with either a <5[Fe/H]/<5rogr ~ —0.3 
dex per decade gradient at a fixed 12.5 Gyr (present) age, 
or with an approximately 3 Gyr per decade age gradient at 
a fixed metallicity of [Fe/H]=+0.09, consistent with most 
previous work on this subject. The age-metallicity degen- 
eracy, of course, prevents us from determining, from these 
data alone, which scenario is most likely. 

We therefore follow previous work (Saglia et al. 2000, 
Tamura & Ohta 2000, La Barbera et al. 2003) in studying 
the evolution of colour gradients as a function of lookback 
time. Our B' , V' and I' bands correspond approximately 
to the rest-frame U, g and R bands. Since few local studies 
are carried out in the g band, we compare our results with 
observations in U — V and V — R. Local values for U — 
V and V — R gradients for cluster galaxies are taken from 
the compilation of Idiart et al. (2002). We use the U — V 
gradients for a z = 0.38 cluster presented in Saglia et al. 
(2000); we have been unable to find V — R data at higher 
redshift. 
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Figure 3. Amplitude of the colour gradients in F450W— F606W 
(top panel) and F606W— F814W (bottom panel) as a function of 
K magnitude. Symbols as in Fig. 2. 



We plot these data and the mean values for colour gra- 
dients in A2218 in Figure 4, together with four models from 
La Barbera et al. (2003) that appear to be most appropri- 
ate to reproduce their observations. The parameters of these 
models are shown in Table 2, where the first column is the 
model ID and the subsequent columns represent the age at 
centre and outskirts, the metal abundance at centre and out- 
skirts and the exponential decay time of the star formation. 
Here, model Tl has a central age of 13 Gyr and an outer age 
of 8 Gyr, with both stellar populations having a metallicity 
of +0.09: it corresponds to the model of the same name in 
La Barbera et al. (2003). Model Zl has two stellar popu- 
larions of the same age (13 Gyr) but a metallicity gradient 
of 0.4 dex from center to edge. Model TZ has both a small 
age gradient (13 Gyr at the center and 11 at the edge) and 
a metal abundance gradient as in model Tl. Finally, model 
T3 corresponds to La Barbera's model ~ Tl (for spirals), 
has a small age gradient (2 Gyr from center to edge, with 
a central age of 13 Gyr) and metal abundance of —0.33 but 
r of 3 Gyr, while all other models have r of 1 Gyr (where 
r is the exponential decay time for star formation). The ac- 
tual values for the ages and metallicities of the models are 
slightly different from those used in La Barbera et al. (2003) , 
as the cosmology and metallicity have changed in the latest 
GALAXEV models. We use the GALAXEV values that best 
approximate those used by La Barbera et al. (2003). 

The models are unable to reproduce the central colours 
(Figure 2), which suggests that these galaxies may be more 
metal rich than [Fe/H]=+0.55 (this is likely due to the typ- 
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Figure 4. Mean colour gradients in B'—V (top panel) and V —V 
(bottom panel) for A2218 galaxies (circles for Es and squares for 
SOs; note the arbitrary offset from the cluster redshift introduced 
for purposes of clarity in presentation). Also plotted are litera- 
ture values (Filled 'up' triangles: Idiart et al. 2002; filled 'down' 
triangles: Saglia et al. 2000) and four models whose parameters 
are described in Table 2 and in the text. 



ical a element overabundances). Because of this, there is 
some doubt as to whether the derived gradients may be 
fairly compared to the observations. The actual values of 
the colour gradients from the models at z = 0.18 are, for 
models T1,Z1,TZ and T3, respectively, -0.40, -0.28, -0.35 
and -0.19 mag. in U~V and -0.08, -0.06, -0.08 and -0.05 
mag. in V — R. While all models are able to reproduce the 
V — R gradients, they are a poor fit to the U — V gradi- 
ent, with the possible exception of model T3. However, if 
we consider the evolution of the gradients, model T3 pre- 
dicts a gradient of —0.13 at z — and —0.30 at z — 0.38, 
which compares badly with the measured values of —0.17 
at z = (Idiart et al. 2002) and -0.08 at z = 0.38 (Saglia 
et al. 2000). Since the models do not appear to reproduce 
the values of the gradients or the galaxy central colours, wc 
choose to carry out a differential comparison by forcing the 
models in Figure 4 to run through the A2218 measurements. 

The main conclusion from this exercise is that the best 
model is the familiar small metallicity gradient (although a 
slightly smaller gradient may perform better), while a small 
age gradient model with large r is marginally worse. Note 
that this latter model performs somewhat better in repro- 
ducing the actual values of the gradients. The combined age 



and metallicity model and the age gradient model are largely 
excluded by the data. 

The similarity in the size of colour gradients for E and 
SO galaxies suggests that their bulges, at least, have sim- 
ilar stellar populations and that these are also similarly 
distributed. The result confirms and extends the observa- 
tions of Ziegler et al. (2001) and Mehlert et al. (2003) for 
bright SO galaxies in Abell 2218 and the Coma cluster re- 
spectively: on the other hand, both Kunstchner & Davies 
(1998) and Mehlert et al. (2003) find that at least some SOs 
in Fornax and Coma, respectively, possess younger stellar 
populations, although these objects tend to be among the 
fainter SOs rather than the bright SOs. Similarly Smail et al. 
(2001) argue that the fainter SOs in their sample (K > 17, or 
L < 0.1L*) may be more distinct from ellipticals. This need 
not be inconsistent with our findings, if bright SOs form in a 
different manner from the faint SOs: in particular, the evo- 
lution may take place mostly in the disks, while our results 
concern mostly bulges. 

Our main conclusions are therefore that the observed 
trends in colour gradients are difficult to explain in the con- 
text of any of the popular models of galaxy formation, but 
stress the essential similarity of the stellar populations of 
bulge dominated systems. Observations of a larger sample 
of objects at higher resolution with ACS or WFPC3 should 
allow us to investigate these issues further. Investigations 
of colour gradients for SO disks should allow us to achieve 
a better understanding of the relation between bulge and 
disk evolution, particularly with reference to the scenario in 
which spiral disks are converted to SO disks via a variety of 
processes (ram stripping, truncation, starbursts induced by 
mergers or interactions, harassment, among others). 
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T 


Tl 


13 


8 
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+0.09 


1 


Zl 


13 


13 
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